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Abstract. We show results from a 7 year VLBI monitoring programme 
of 3C 273 at millimeter wavelengths. We find evidence for component 
■ acceleration, motion or rotation of fluid dynamical patterns, and an 

outburst-ejection relation between 7— ray flares and new jet components. 
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^r^ ' 1. Introduction 



VLBI imaging at millimeter wavelengths (mm- VLBI) gives highest angular res- 
olution and facilitates direct studies of the self-absorbed jet base in AGN. For 
(3 I 3C273 {z = 0.158) images with an angular resolution of up to 50 //as (1 //as 

= 10~^ arcsec) are made at 86 GHz. This corresponds to a spatial scale of 
~ 1000 Schwarzschild radii (for M = 10^ Mq ). The small observing beam 
• allows to accurately determine positions of jet components. This facilitates de- 

C^l tailed studies of the jet structure and kinematics near the nucleus, in particular 

Q ■ with regard to the broad-band (radio to Gamma-ray) flux density variability 

^ I and the birth of new 'VLBI components'. Here we summarize new results from 

^2 ■ our multi epoch (1990 - 1997) and high observing frequency (15, 22, 43 and 

86 GHz) VLBI monitoring. 



2. Results 

. At sub-milliarcsecond resolution, 3C 273 shows a one sided core-jet structure 

of several milliarcseconds length. The jet breaks up into multiple VLBI com- 
ponents, which - when represented by Gaussian components - seem to sepa- 
rate at apparent superluminal speeds from the stationary assumed VLBI core. 
The cross-identification of the model-fit components, seen at different times and 
epochs, is facilitated by small (< 0.2 mas), and to first order negligible, opacity 
shifts of the component positions relative to the VLBI-core. Quasi-simultaneous 
data sets (cf. Fig. 1) demonstrate convincingly the reliability of the compo- 
nent identification, which results in a kinematic scenario, in which all detected 
jet components (C6 - C18) move steadily (without 'jumps' in position) away 
from the core (Fig. 4). For the components with enough data points at small 
(< 2 mas) and large (> 2 mas) core separations, quadratic fits to the radial mo- 
tion r(t) (but also for x(t) -right ascension, and y(t) -declination) represent the 
observations much better than linear fits. Thus the components seem to accel- 
erate as they move out. The velocities range typically from Papp = 4 — 8 (for 
Ho = 100 km s"^ Mpc-\ qo = 0.5). 
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Fig. 1: 3C273 at 22 GHz (top), 43 GHz (cen- 
ter), and 86 GHz (bottom) observed in 1995.15 
- 1995.18. Contour levels are -0.5, 0.5, 1, 2, 5, 
10, 15, 30, 50, 70, and 90% of the peak of 3.0 
(top), 5.4 (center), and 4.7 Jy/beam (bottom). 
For the 22 GHz map, the 0.5 % contour is omit- 
ted. The restoring beam is 0.4 x 0.15 mas in 
size, oriented at pa = ° . The maps are cen- 
tered on the eastern component (the core), the 
dashed lines guide the eye and help to identify 
corresponding jet components. 
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Fig. 2: Spectral index variations along the jet. 
For 1997 (circles, solid lines), the spectral in- 
dex gradient is calculated directly from the in- 
tensity profiles of the maps at 15 and 86 GHz. 
For 1995 (squares, dashed line), the spectral 
indices were derived from Gaussian component 
model fits at 22 and 86 GHz. We note that 
during the time interval 1995 - 1997, different 
jet components occupied this jet region. The 
spectral profile along the jet, however, did not 
change significantly. 
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Fig. 3: The motion of the mean jet axis in the 
inner jet of 3C 273 at 15 GHz. Symbols denote 
for different epochs: 1995.54 (circles), 1996.94 
(squares), 1997.04 (diamonds), and 1997.19 
(triangles). The transverse oscillations of the 
ridge line are measured relative to a straight 
line oriented along pa — 240 ° . We note the 
systematic longitudinal displacement of max- 
ima and minima. This corresponds to an ap- 



parent pattern velocity of 0a 



: 4.2, which is 



by a factor of up to 2 slower than the compo- 
nent motion. 
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Dual-frequency maps obtained in 1995 (22/86 GHz) and 1997 (15/86 GHz) 
allow to measure the spectral index gradient along the jet (see Fig. 2). The 
spectrum oscillates between —1.0 < a < +0.5 {Si^ oc z/°). Most noteworthy, the 
spectral gradients did not change significantly over the 2 year time period, al- 
though different jet components occupied this jet region. Thus, the geometrical 
(eg. relativistic aberration) and/or the physical environment (pressure, density, 
B-field) in the jet must determine the observed properties of the VLBI compo- 
nents. Hence, the latter do not form 'physical entities', but seem to react to the 
physical conditions of the jet fluid. 

Further evidence for a fluid dynamical interpretation comes from a study 
of the mean jet axis and the transverse width of the jet. Both oscillate quasi- 
sinusoidally on mas-scales. At 15 GHz, the variation of the ridge-line with time 
could be determined from 4 VLBA maps obtained during 1995 - 1997 (Fig. 3). 
The maxima and minima of the ridge-line are systematically displaced. This 
'longitudinal' shift suggests motion with a pattern velocity of (3app = 4.2. The 
sinusoidal curvature of the jet axis, however, is more indicative for jet rota- 
tion rather than for longitudinal waves. Helical Kclvin-Hclmlioltz instabilities 
propagating in the jet sheath could mimic such rotation, which, when seen in 
projection, would explain also the spectral index oscillation in Figure 2. 

For many AGN, a correlation between flux density variability and ejection 
of VLBI components is suggested. In 3C 273 we identified 13 jet components 
(C6 - C18) and traced their motion back to their ejection from the VLBI core. 
The typical measurement uncertainty for the ejection times to ranges between 
0.2 - O.Syr. In Figure 5 (left) we plot to and the millimeter-variability (22 - 
230 GHz). We also add the Gamma-ray detections of 3C 273 from EGRET. In 
Figure 5 (right) we plot the onset times of the mm-flares derived from these 
light curves (Tiirler et al. 1999) together with the VLBI ejection time {to) and 
the Gamma-ray fluxes. For each onset of a mm-flare, we find that a new jet 
component was ejected. 
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Although the time samphng of the Gamma-ray data is quite coarse, a 
relation between component ejection and high Gamma-ray flux appears very 
likely (note that each Gamma-detection already means higher than usual 7- 
brightness). From a more detailed analysis (Krichbaum et al. 2001) we ob- 
tain for the time lag between component ejection and onset of a mm-flare: 
to — t™™ = 0.1 lb 0.2 yr. If we assume that the observed peaks in the Gamma- 
ray light-curve are located near the times of flux density maxima, wc obtain 
^0 ~^o™ = 0.3 lb 0.3 yr. Although the Gamma-ray variability may be faster, this 
result is fully consistent with the more general finding of enhanced Gamma-ray 
fluxes mainly during the rising phase of millimeter flares. We therefore suggest 
the following tentative sequence of events: t™™ 1^ to 1^ - the onset of a mil- 
limeter flare is followed by the ejection of a new VLBI component and, either 
simultaneously or slightly time-delayed, an increase of the Gamma-ray flux. If 
we focus only on those VLBI components, which were ejected close to the main 
maxima of the Gamma-ray light-curve in Figure 5, we obtain time lags of tg — 
of < 0.5 yr for C12, < 0.9yr for CIS, < 0.2 yr for C16 and < 0.1 yr for C18. In 
all cases the Gamma-rays seem to peak a little later than the time of component 
ejection. With Papp — 4 near the core, the Gamma-rays would then escape at 
a radius < 0.1 mas. This corresponds to < 2000 Schwarzschild radii (for 
a 10^ Mq black hole) or < 6 • 10^''' cm, consistent with theoretical expectations, 
in which Gamma-rays escape the horizon of photon-photon pair production at 
separations of a few hundred to a few thousand Schwarzschild radii. 
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Fig. 5: Left: Flux density variations at 230 GHz (filled diamonds), 86 GHz (open circles) 
and 22 GHz (filled squares). Upward oriented triangles denote Gammarray fluxes from 
EGRET, downward oriented triangles arc upper limits. The extrapolated ejection times of 
the VLBI components and their uncertainties are indicated by filled circles with horizontal 
bars along the time axis. Right: Broad-band flux density activity and component ejection. 
VLBI component ejection (open circles), Gamma-ray fluxes (triangles, downward oriented for 
upper limits) and onset-times for the millimeter flares (open squares, from Tiirler et al. 1999). 
Labels denote the component identification from VLBI. 
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